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ABSTRACT 


The principal features of this document are the analysis of a large dual-reflector 
antenna system by vector Kirchhoff theory, the evaluation of subreflector 
aperture-blocking, determination of the diffraction and blockage effects of a 
subreflector mounting structure, and an estimate of strut-blockage effects. 

Most of the computations are for a frequency of 15.3 GHz, and were carried out 
using the IBM 360/91 and 360/95 systems at Goddard Space Flight Center. 

The Fortran IV computer program used to perform the computations is of a 
general and modular type so that various system parameters such as frequency, 
eccentricity, diameter, focal-length, etc. can be varied at will. The parameters 
of the 60-foot NRL Ku-band installation at Waldorf, Maryland, were entered into 
the program for purposes of this report. Similar calculations could be performed 
for the NELC installation at La Posta, California, the NASA Wallops Station 
facility in Virginia, and other antenna systems, by a simple change in IBM con- 
trol cards. A comparison is made between secondary radiation patterns of the 
NRL antenna measured by DOD Satellite and those obtained by analytical/numeri- 
cal methods at a frequency of 7.3 GHz, 
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GLOSSARY OF NOTATION 


Symbol 

Meaning 

k 

wavelength 

cpu 

central processor unit 

ATS 

Applications Technology Satellite 

GHz 

giga-Hertz 

E(x', y', z'), H(x', y\ z') 

scattered electric and magnetic vector fields 

x, y, z 

Cartesian coordinates 

r, 9, 4> 

spherical coordinates 

fi i 

unit normal to surface 

0 

Solution to wave equation 

V \fj 

vector gradient of 

ds 

differential area 

Mi 6 

constitutive parameters: magnetic permeability; 
inductive capacity 

O) 

angular frequency 

j 

imaginary operator = 


electric and magnetic polarization vectors 

k 

wave number 

AAA 

f ^ J 

spherical basis vectors 

A A A 

i» j* k 

Cartesian basis vectors 

F 

focal length 


surface 


radial and angle variables of a cylindrical net 

s 

source strength 

N 

exponent of directive source 

5 

directive point source 

V 

polarization moment 
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Symbol 


Meaning 


a-,/3, y 

Euler angles 

P e 

feed displacement vector with components (x € , y e , z f 

D 

diameter 

a, c 

parameters of a hyperboloid 

LIj, LI 2 

sampling resolution on surfaces y t , y 2 respectively 
in wavelengths 

EOS 

Earth Observation Satellite 

R, 0 , $ 

spherical coordinates for observing feeds and sub- 
systems 

A 

area 

z lf z 2 

complex numbers 

BW 

beam width 

SLL 

side-lobe level 

db 

decibel 


reference plane 

A 

increment 

f 

frequency 

4> 

phase (note context) 
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DIFFRACTION STUDIES APPLICABLE TO 
60-FOOT. MICROWAVE RESEARCH FACILITIES 


INTRODUCTION 

The Antenna Systems Branch, Network Engineering Division at Goddard Space 
Flight Center has developed, and maintains, a general-purpose digital computer 
program for calculating diffraction patterns of reflector- type antenna systems. 

A dual- reflector capability was achieved about one year ago to compute both 
Fresnel and Fraunhofer diffraction patterns in either the transmit or receive 
mode of operation. The program integrates over charge and current distribu- 
tions on reflector surfaces in ordinary physical 3 -space to produce scattered 
electric and magnetic field magnitudes, wavefronts, and time-average Poynting 
vectors, using a complex -vector Kirchhoff approach with Kottler boundary in- 
tegral, satisfying Maxwell's equations. 

Discussions between GSFC, NRL, and NELC personnel indicated that the subject 
diffraction program was well-suited for providing the following characteristics: 
directive gain, aperture efficiency, spillover efficiency, diffraction effects from 
mounting plates, and an estimate of secondary blockage effects due to the hyper- 
boloidal subreflector and the associated support struts for large microwave 
facilities at K-band. Data on beamwidth, side-lobe levels and cross-polarization 
components are also provided by the program. A point of concern at the outset 
of this effort, which included a comparison of computed results with actual 
measured patterns obtained via DOD Satellite at 7.3 GHz, was the electrical size 
of the main aperture at 15.3 GHz. The GSFC program had given excellent results 
for dual-reflector systems with reflectors around 100 \ by 100 and single- 
reflector systems with reflectors approximately 500 X by 500 \, but the present 
study of the Waldorf antenna, Figure 1, involved an aperture approaching 1000 \ 
by 1000 A. and was also a dual- reflector configuration.* 

Computation of the diffraction patterns of large K-band dual reflector antennas, 
particularly when appearing in array, had presented a challenge for some time 
in TDRSS studies. The analysis of the 60-foot NASA/NRL/NELC microwave in- 
stallations was, therefore, undertaken to determine whether or not it was feasible 
to obtain meaningful results for such a problem with the existing integration 
algorithm. Annexation of the strut-blockage subroutine via a Babinet’s principle 
was an added complication that had received only limited verification at the out- 
set of this study. The composite subreflector posed no difficulty as this capability 


*The NRL Waldorf installation was funded by NASA to collect ATS- 5 rainstorm data at 15 GHz. 
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had been tested extensively several years earlier. Some modification of the 
basic program was anticipated — and required — due to the inevitable trade-off 
between computer cpu-time and accuracy. Computer program development is 
usually determined by actual requirements and applications, and in this instance 
direct benefits were realized for EOS 3-reflector passive multichannel micro- 
wave radiometer (PMMR) antenna computations, ATS-F diffraction studies, 
STDN 40-foot and 85-foot antenna analyses, and TDRSS trade-off studies for 
determination of the antenna payload. 


FORMULATION OF THE SCATTERED FIELDS 

hi the transient mode of operation the fields due to the main or parabolic reflector 
( 7 2 ) of the system are obtained from 

E(x, y', z') = - — j* [(n xH 2 )- V] ds - jco/i. ^ J (n 2 xH 2 ) 0 ds, 

7 2 y 2 


H(x\ y', z) 


1 

477 



(n 2 x H 2 ) x V \p ds, 


where 


0 


e-i>” 


r 


and 


= - ( jk + 1 )+ i r , 


taking l r in a local context on the reflecting surface. 

Since only the magnetic field H 2 is required on y 2 , the integral 

H 2 (x 2 , y ' 2 , z'i) = - 4 — J (n x xHj) x Vv/» ds 

• y \ 
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over the hyperbolic subreflector (y x ) is sufficient to, obtain the illumination of y 2 
due to y l . The subsystem surface maybe extended to include the flat mounting 
plate {y\). 

The preceding formulation is applicable in both far-field (Fraunhofer) regions 
and intermediate near-field (Fresnel) regions, and is well-suited tp dual or 
multiple-reflector systems, where a point on one surface may lie in the Fresnel 
region of another surface over which an illumination distribution is being inte- 
grated. It can be shown that the formulation used here is identical to the 
Kirchhoff-Kottler and Franz formulations, but is better-suited to numerical 
integration methods than the latter.* 


GEOMETRY 

The geometry for the problem is shown in a photograph designated Figure 1. 
Input parameters and surface equations are as follows. 

Paraboloid (y 3 ) 


x - cr sin £ , y - - a cos £ , z 



z ip 


z,p = -F = -18.0’ 
a = 30.0’ 

max 

cr Q = 0.0' (initially, without obscuration) 

Hyperboloid (y t ) 

x - cr sin £ , y = - a cos £, z ~ c ( 1 + cr 2 / a 2 ) + z 1H 

Z 1H = -ce = -5.87’ 
a = 4.03’ 
c = 4.27’ 

*See Ref. 1, 2, 3, and 4. 
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3.00' 


C r 

max 

= o.o- 

Mounting Plate (rl) 

x = <y sin £, y ~ -a cos £, z = c c cr + z lc 

z lc = 0.0' 

C c = 00 

= 3,00’ 
cr = 4 41* 

Feed (?) 

? = Scos N ® =*■ directive field 
S = 1.0 =» source strength 

N = 33.2 approximately -10 db edge taper on y x 
P(I) = 1.0, P(J) = P(K) = 0.0 =► polarization moments 

a = /3 = y = 0,0 =*• Euler angles for feed 

x e - y £ ~ 0.0’, z e ® -2c e = -11.7' => feed location 

F/D = 0.30 =>• F - number of paraboloid 


The subroutines for the paraboloid and the hyperboloid require no elaboration, 
however, the subroutine for the flat mounting plate deserves a word of explana- 
tion. hi the physical installation the plate is flat and is situated in the focal- 
plane of the paraboloid. From the point of view of the simulation, the electro- 
magnetic scattering from the plate is due to that portion of the plate which is not 
shielded from the feed by the hyperboloid. The plate can therefore be thought of 
as having a central hole of radius a- = 3.00’, the radius of the hyperboloid. 

Although the plate is square, 93.46 inches along an edge, the equivalent area of 
a flat circular disc was used to compute the effects of the mounting plate. To be 
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more precise, the annulus of a flat circular disc was used in the simulation. 

This permitted the utilization of an existing subroutine (cone) which exhibited 
rotational symmetry. All integration over y 1 and y[ was then done by means 
of a cylindrical (cr, £) net, and a Cartesian (x, y) net for the square plate was 
avoided. The parameter c c is the inverse-slope of a cone; the latter degenerates 
into a flat disc when c c vanishes. Choosing z lc - 0.0' situates the flat disc in 
the focal plane. 


RESULTS OF THE "LITERAL” METHOD 

Under ray optics, an ordinary Cassegrain system constitutes a one-to-one map- 
ping from the sub reflector onto main reflector: Under diffraction, the illumina- 
tion of the main reflector is obtained from a many-to-one mapping of the illumi- 
nation distribution (charge and current) of the subreflector onto the main 
reflector. This many-to-one mapping has been termed the "literal” method 
since it is the natural result of redundantly applying the Huygens principle to 
multiple-reflector systems. The reader may wish to assess the awkwardness 
of this approach when three or more surfaces make up the system.* 

i 

Since application of the "literal" method was expected to result in excessive 
cpu-time and core requirements for a dual -reflector system in which the main 
reflector exceeded 800 k by 800 k , an approach which regarded the prime-feed/ 
subreflector combination as an effective "equivalent source" was favored. Fig- 
ure 2 shows the far-field pattern via the "literal" method for this Cassagrain 
system without mounting-plate, central obscuration, or strut obscuration at a 
frequency of 1.53 GHz. Results appear acceptable out to 3 degrees in the theta 
domain for this low-frequency test case obtained at a sampling of LIj = 0.25 
and LI 2 = 2.5 on the subreflector and main reflector respectively. The far-field 
pattern at 15.3 GHz was also obtained via the "literal" method, using the same 
input parameters as before, and is shown as Figure 3. This result was obtained 
for LIj = 1.25 and LI 2 =? 6.00, storing the illumination values on y i and y 2 , and 
taking advantage of the biplanar symmetry that exists for the problem with re- 
spect to illumination on the reflector surfaces. Results appear credible out to 
0.3 degrees in the theta domain although null-filling is in evidence at N 4 . (Null- 
filling is not certain at N 2 since Ad - 0.005 degree does not resolve a good 
minimum or true null unambiguously.) Considering that a ten-fold increase in 
frequency occurred between Figure 2 and Figure 3, the side-lobe structure out to 
0.3 degrees. corresponds to ■■3.0. degrees of the previous computation. The onset 


*A three-reflector "telescope" system was studied for the Earth Observation Satellite (EOS), 
however, the ''literal” method was avoided due to prohibitive cpu-times involved. 
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of failure of the solutions due to inadequate integration sampling (Lip LI 2 ) is 
apparent from the null-filling, anomalous null-to-null sidelobe angular width, 
sidelobe level, and departure of side-lobe peaks from an asymptote or monotonic 
law of descent. It can be seen from Figure 2 and Figure 3 that degradation sets 
in progressively for larger values of theta ( 6 ). 

Since the practical limits of cpu time and core had already been exceeded in 
obtaining results at 15.3 GHz, Figure 3, and since the effects of the flat mounting 
plate and struts had not yet been included, the "literal" method was abandoned 
for the remainder of this particular study. The "equivalent source" method 
which derives scattered field data from the prime-feed subreflector combination 
was adopted, and even this approach was found to be demanding with respect to 
cpu-time and core when the effects of the mounting plate were included. The 
annular area of the mounting plate was found to be about equal to that of the 
hyperboloid. Integration sampling intervals somewhat greater than those pre- 
scribed by the antenna analogue to Shannon's sampling theorem* were finally 
used, but this precedes the development. 


SUBSYSTEM PATTERNS 

The "equivalent source" method is now detailed without inclusion of the annular 
mounting plate to illustrate some of the problems and show the reader the basis 
for the sampling (LI X ) on the subreflector (y x ). First of all, there is the question 
of the range (R) at which field values scattered by the prime^feed/ subreflector 
combination should be reevaluated. See the expression for H 2 (x 2 , y 2 , z 2 ). 
Clearly, the values on a parabolic locus (y 2 ) are ultimately required, but this 
view obscures the fact that the prime-feed/hyperboloid subsystem behaves as a 
virtual spherical point- source at F. A large value of R, to simulate infinite 
range, has some merit. It is possible, however, that in general there may be 
near-field considerations since any point on y 2 , at 15.3 GHz, is in the interme- 
diate near-field of if the physical dimensions of the latter are regarded. 
Another viewpoint is that the effective diameter for a hyperboloid fed at the con- 
jugate focal point is not anything like the physical diameter, but is essentially 
zero (i.e. a point source) in R = 2d 2 A, the usual criterion for discerning "near- 
field" and "far-field". 

The range (R) was initially chosen equal to the focal distance (F) of the paraboloid 
for the "equivalent source" method since this allows one to observe the purity of 
the wavefront on a circular arc whose center is the phase center of the prime- 
feed/hyperboloid subsystem being studied. This is especially convenient because 


*Ref. 5, page 321. 
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the nearly constant phase on this spherical arc (R = F) can now be used to 
further develop the "equivalent source" method. Both field phase and amplitude 
on this arc are placed in core for a preselected resolution in theta (A 6 = 2°) 
between the vertex of y 2 and the geometrical bound or edge of y 2 . This "equiva- 
lent source" data is then modified by the factor e _ i kp /p which effectively pro- 
vides for the differential phase between the vertex wavefront (R = F) and the 
paraboloid ( y 2 ) as well as the space divergence of the amplitude of the illumina- 
tion arriving at y 2 . This is in lieu of computing the difference between the arc 
R = F and p to the paraboloid for each differential (incremental) area. 

The absolute phase of the source has then been changed in the amount e _jkF and 
the original source strength has been multiplied by 1/ F, both of which are 
irrelevant with regard to the final diffraction pattern. This can be seen as follows. 

Assume a virtual point source of strength S 0 , initial phase 0 O , and directivity 5 
such that 




e 


-jkF 

F 



spherical wave expanded 
out to y 2 


spherical wave expanded 
out to R = F 


where T is the transformation which rigorously preserves both field phase and 
amplitude in going from the arc R = F to the reflector y 2 . 



Now the first factor of T was used, as stated previously, but the second factor 
was ignored (i.e., read-in "unity"). Failure to include the second factor there- 
fore multiplies the assumed source by the reciprocal of that factor; 
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so that the omission results in nothing more than a reinterpretation of the real 
source as 



j (0O-kF) cf 


In other words, amplitude and initial phase of the assumed source have been 
modified in an irrelevant manner. The validity of the preceding discussion is 
further borne out by a comparison of the secondary radiation patterns by the 
"literal” method with those of the "equivalent-source" method using the prime- 
feed/hyperboloid subsystem. Obviously, the approach is predicated on a spheri- 
cal law of divergence between the arc R = F and the value p on the main re- 
flector (y 2 ). 

Figure 4 shows the (relative) phase and amplitude of the totality of the Huygens 
wavelets scattered from the hyperboloid (7 1 ) onto the spherical arc (R = F = 18') 
passing through the vertex of y 2 , for a sampling interval of LI t = 1.25, Note the 
geometrical bound of y 2 . The result is extremely coarse, however, .surprisingly 
good secondary patterns can sometimes be obtained from the main reflector 
with such an "equivalent source". Phase dominates problems of this type, and 
rather large amplitude variations are a lesser consideration. In Figure 4 the 
phase excursions approach ±180 degrees, and it is remarkable that good results 
can be obtained for such departure from a pure spherical wavefront. 

Figure 5 shows the (relative) phase and amplitude of the scattered electric field 
on R * F = 18’ as for Figure 4, but LI t = 0.600. It can be seen that the large 
"spike" in amplitude has diminished near 0 = 180°, and the variations, generally, 
are smaller. More importantly, the phase has become constant to within ±20 de- 
grees or so for 100° < 0 5 140°, which represents considerable energy in an 
annular region of the sphere R — F = 18'. The phase variations still approach 
±180 degrees in the domain 140° < 0 < 180°. Also note the behavior of the phase 
plot outside of the geometrical bound 90° < 0 < 100°. It is still violently oscil- 
latory, but less so than in the preceding instance. 

Figure 6 shows the (relative) phase and amplitude of the scattered field in R = 

F = 18' for LI t = 0.300, the sampling used in most subsequent dual -reflector 
computations. The improvement over Figure 5 is striking. Both phase and 
amplitude closely approach the classical ray optics concept for scattering by the 
subsystem. Over the domain of interest, 100° 5 ® 5 180°, the amplitude is a 
monotonically decreasing function with a ripple of approximately ±1 db, while 
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the phase plot or wavefront exhibits a variation of about ±10 degrees. The scat- 
tered values of Figure 6 were considered "stable" solutions with respect to the 
integration sampling (LI X ). Once stored, these values are not computed again 
even though main reflector obscuration by the hyperboloid and/or the struts is 
considered subsequently. None of these values are computed anew for different 

secondary pattern cuts 4> = 0°, <p = 45°, 4> = 90% etc. 

\ 

Figure 7 shows the subsystem phase and amplitude pattern at K = F = 18’ with 
the flat annular mounting plate included. The reader will notice that the hyper- 
boloidal part of the subsystem was sampled at LI X = 0.300, but the flat annular 
mounting plate was sampled at LI ' = 1.00 in carrying out the integration process. 
It was determined by means of separate calculations with the program that the 
differences between LI x r = 0.500 and LI/ = 1.00 were negligible, therefore, the 
coarser interval was selected in forming this subsystem pattern in the interest 
of cpu-time reduction. 

The phase and amplitude of the backscattered fields in Figure 7 are perturbed 
by the influence of the flat annular disc.* In the region 100° < 0 < 160°, which 
is between the geometrical bounds of y x and y 2 , the amplitude exhibits a normal 
monotonic descent over most of the area of interest. The exception is the interval 
140° 166°. The phase is constant, as required, to within ±15 electrical 

degrees or so over most of the region of interest, .however, it varies ±50 electri- 
cal degrees in the excepted interval. Confining the discussion to the interval 
100° < 0 £ 166° between the geometrical bounds, and regarding areas on a unit 
sphere through which power is flowing, the percentage of perturbed to unperturbed 
power is estimated (without respect to power density variations) by areas alone. 


11 - ■ 11 


A = lids = I I R 2 sin 0 d 0 d = - 2 tt cos 0 1 2 

' 0 , 


At 0 tal " “277 (cos 166° - cos 100°) 




ERTURBED 


- - 2v (cos 166° - cos 140°) 


\oTAL 




X 100 = • .25% 


ERTURBED 


*See Appendix B. 
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It is anticipated, therefore, that the Kirchhoff integration over the illumination 
distribution on y 2 > using Figure 7 data, will not greatly affect the far-field radia- 
tion pattern. 

The polarization of the "equivalent source" enters the problem when the integral 
H 2 (x 2 ', y 2 ' , z j ) is evaluated at R = F. It should be noted that the field quantities 
are complex— vectors due to the unit polarization vectors (hj) of the real or 
physical source. Let the components of the source polarization moment (desig- 
nated P(I), P(J), P (K) in the program input) be (v x , v y , v z ) = v. Then the mag- 
netic polarization vector incident or 7j is 

- _ p_ X V 

hi Ip x V I 


and the electric polarization vector is 


P X (y up') 
Ip X (v xp)| 


Only hi is required in the present computations, although e x is probably more 
familiar and more easily visualized, since only incident magnetic fields H t = 
|H, |h t appear in the diffraction integral 


“2 C x 2’ y 2’ Z 2 ) 



X H 1 ) X V0 ds 


at R = F. 

The "equivalent source" method is relatively straightforward, with one possible 
exception. It has been shown that the backscattered illumination is a set of 
complex-vectors, and the vector character of the problem enters via the unit 
polarization vector (hj). The complex character of the problem is associated 
with the time domain or phase of the fields at various points. Although the present 
computations were not affected, a delicate question arises concerning intermedi- 
ate complex values for analogue data. For example, the E -plane, H-plane eccen- 
tricity of the backscattered fields were disregarded here, and the spherical cut 
(r, 6, 4>) = {18 \6 , 0°) provided phase and amplitude data for all other cuts via a 
trivial interpolation between 4> - 0° and 0 = 360°. Suppose the pattern 
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eccentricity had been large. Additional cuts at = 45° or A 0 = 15°, etc. 
could have been employed to display the character of the backscattered fields. 
This would have raised the following question. 

What is the "average” of two complex numbers (ordered pairs)? How does one 
interpolate or find an intermediate complex number between two given complex 
values? One approach might be to take the average of the magnitudes and the 
average of the arguments or phase angles. Another approach might .be to aver- 
age the real and imaginary parts separately to produce an intermediate value. 
Unfortunately, the two approaches lead to different results for an arbitrary 
choice of "nearly equal" complex numbers: z lt z 2 ■ Both methods can be very 
disappointing from a practical point of view, not to mention the rigorous point of 
view. It is possible to obtain such interesting results as an argument or phase 
which is approximately 77 radians different from the argument of z l and z 2 . 

Also, it is possible to obtain an amplitude which is either literally zero or ap- 
proximately zero even though the magnitudes of z 1 , z 2 were arbitrarily large I 
An area of difficulty has been identified here, however, the impact on the present 
problem is negligible. Many "cuts" on analogue data relieve the severity of the 
problem as do a nearly spherical wavefront and small E-plane, H-plane eccen- 
tricity. 


RESULTS OF THE "EQUIVALENT SOURCE" METHOD 

The reference pattern , in the following sequence of 15.3 GHz far-field radiation 
patterns from the dual-reflector system, is Figure 8 . As stated previously, sub- 
reflector sampling was stable at LIj = 0.300. It was found that main reflector 
sampling was stable at LI 2 = 6.00 if the first three or four sidelobes only were 
required for study.* Beamwidth <BW) equals approximately .070 degree, and the 
first sidelobe level (SLL t ) is approximately -25.5 db below the beam maximum. 

The first departure from the reference pattern is Figure 9, which shows the 
effect of central obscuration due to the hyperboloid ( a 0 = 0.0', a max = 3.0'). 

The range of integration over y 2 is now restricted by cr 0 = 3.0', a max = 30.0'. 

A weak interferometer effect is anticipated and observed. (BW ~ .070 degree, 
SLLj ~ -23.0 db.) A discussion of directive gain degradation is deferred to a 
summary at the end of this section. 


*Grating lobes of classical array theory do not appear due to the spatial distribution of the LI 
samples on the polar (a, £) net. 
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The second departure from the reference pattern is Figure 10, and shows the 
effect of central obscuration due to the hyperboloid and the flat mounting plate 
bounded by cr Q = 0.0', cr max = 4.41'. The range of integration over y 2 now 
restricted by cr 0 = 4.41 1 , o- max = 30,0'. In this calculation the backscattering 
from the annular mounting plate is also included. Recall Figure 7. It can be 
seen that beamwidth is negligibly affected relative to the preceding Figure 9, 
but the sidelobe structure has been altered significantly. (SLL x % - 19.5 db.) 

The third departure from the reference pattern is Figure 11, and is the "com- 
plete" simulation of the physical structure for this report. Hyperboloid and 
mounting plate obscuration as well as scattering are present as for Figure 10, 
but strut blockage has now been taken into account. Beamwidth is not altered 
significantly, and has narrowed slightly, if anything. Sidelobe structure has 
taken on some unusual characteristics. Pattern level changes are not very large. 
(SLLj % - 19.2 db.) The 4> = 90° cut,\not shown, closely resembled the 4> = 0° 
cut of Figure 11 in all except the fine detail. A 4> = 45° cut, also not shown, ex- 
hibited relatively high sidelobes (SLLj % - 17.5 db). 


The level of the electric field E e (main polarization) for the preceding calcula- 
tions, Figure 8 through 11, as derived from the Kirchhoff integration compare 
as follows. 


Figure 8 reference pattern 

Figure 9 hyperboloid obscuration included 

Figure 10 hyperboloid and mounting-plate 
obscuring and scattering 


Figure 11 complete system with (secondary) 
strut blockage 


-94.529 db 
-94.683 db 

- 94.726 db 

- 95.030 db 


The implication of these values is that elements of the system degrade the gain 
approximately as shown below. 


hyperboloid obscuration 
mounting-plate effects 
spar obscuration 


0.15 db 
0.04 db 
0.30 db 


It appears that the combined gain degradation is about 0.5 db and the mounting 
plate, though comparable in area with the hyperboloid, has a negligible influ- 
ence. Details of the development leading to the obscuration estimates can be 
found in the Appendix A of this report. Gross polarization components are a nat- 
ural by-product with the GSFC diffraction program since the vector integration 
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in Cartesian coordinates leads to fields which are simply projected onto the 
spherical basis vectors to produce main and cross-polarization components at 
the close of the calculations. None of the cross -polarization data are presented 
here as no degradation problem was noted in this characteristic of the system. 

The preceding results demonstrated that the existing GSFC diffraction program 
was capable of dealing with large dual-reflector Ku-band antenna systems. A 
credible beam solution was obtained for the 60-foot Cassegrain installation and 
such considerations as cpu-time, core-requirements, significant digits, accuracy, 
etc. did not individually or collectively present any insurmountable problems. 

In a sense, the body of this report represented a first milestone. A secondary 
objective was the correlation of the computed results with a good Ku-band field 
measurement for precisely the same system that had been simulated. This ob- 
jective was not realized since reliable 15.3 GHz data for these sites did not be- 
come available, even though considerable interest existed concerning the shape 
of a Ku-band pattern that was influenced by subsystem obscurations, a large 
mounting-plate, and strut blockage. Measured field-pattern data did become 
available at a lower frequency, however, and it was decided to simulate this con- 
figuration also by means of the diffraction program to obtain at least one com- 
parison between the analytical/numerical approach and field practice. 

The low-frequency (C-band) comparison is detailed in Appendix B of this report. 

It was determined that certain deficiencies existed in the initial analysis when 
the mounting-plate was included, but only the far-out sidelobes were significantly 
affected. Appendix B also includes the results of a 15.3 GHz simulation using the 
improved analysis. 


SUMMARY 

This report has attempted to illustrate some of the practical problems and solu- 
tions for a dual- reflector geometry that was complicated by a compound sub- 
reflector, a central aperture obscuration, and strut blockage of the main aperture 
Perhaps the outstanding obstacle was that the paraboloidal reflector dimensions 
for this Cassegrain system, in wavelengths, precluded utilization of the "literal" 
method of Huygens under a Kirchhoff-Kottler type of formulation. A fundamental 
difficulty concerning the interpolation of complex numbers was encountered. In 
addition, it was shown that expansion of wavefronts from a spherical locus cannot 
proceed by means of a single, simple\spherical divergence rule when multiple 
phase-centers are encountered!.. The degradation of results in such cases will 
depend on the location and weight of those phase centers. An estimate of strut 
blockage effects was included, and this was based on a simple Babinet principle 
which deleted umbral regions for the Kirchhoff integration. It was noted that 
shadow updating is not required when the sidelobe structure of interest lies in a 
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small angular region away from the system axis. Primary blockage, as viewed 
from the virtual origins of the subsystem, was not considered. Mutual coupling 
effects were, likewise, not considered. 

The subjects of system gain and directive gain were not covered comprehensively 
since the principal objective was the computation of pattern shape or form. Ob- 
viously the application of the definition of directive gain 

G = F ( 0 , 0 ) 

4 ^ JJ P ( 0 , 4 >) sin 5 d£ d 4 > 


over the half-space containing the radiation pattern of this system would be ex- 
tremely costly in cpu-time. A comparison method, based on the notion of a de- 
generate paraboloid (F = °°) of identical diameter and uniform phase, amplitude, 
and polarization illumination might be employed successfully here. An estimate 
of gain degradation due to hyperboloid obscuration, mount ing-plate effects, and 
span obscuration was included. A linear polarization state was assumed through- 
out for convenience . 

The 7.3 GHz satellite data for the Waldorf, Maryland installation, as furnished by 
the U.S. Naval Research Laboratory, greatly stimulated interest in the simula- 
tion. A valuable comparison between "theory” and "practice" lead to the explora- 
tion of several methods for computing the \ diffraction patterns. All approaches 
used here simulated the main beam and first sidelobe level fairly well. Signifi- 
cant differences between methods lay in the far-out sidelobe region. Techniques 
evolved in this study are being applied to diffraction problems for ATS-F, STDN, 
EOS, and TDRSS. 

Diffraction studies are being continued, however, the investigation of the present 
configuration is terminated. This is an interim report. 
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appendix a 


SECONDARY BLOCKAGE CALCULATIONS 


Figure 12 depicts a general secondary blockage or obscuration example 0 2 . 

For high-gain, small beamwidth problems the interesting features of the radia- 
tion pattern do not lie far from 6 = 0° in the (9-domain, and the square blockage 
area can be taken as without serious error. An accurate projection of the 
support trusses is not attempted in Figure 12 which illustrates only the general 
nature of the problem and its solution. The obscuration Oj can be taken to 
represent the hyperboloid and square mounting plate encountered earlier in this 
report. In evaluating the Kirchhoff-Kottler class diffraction integral over y 2 it 
is sufficient, for a first-order estimate of blockage effects based on a Babinet's 
principle, to delete the area O t in the integration for the scattered field. 


E (x', y\ z') - 


1 

3 CO £ 


4 7 T 1 

Jv - n 


H 3 ) ■ V V0ds-jo>|U 


^ f 

J'y - 


(n 2 x H 2 ) ds 


In cases where lower-gain, larger-beamwidth problems are encountered, the 
same Babinet's approach requires an up-dating of the shadow for every observer 
position (»= , 9, <p) - (R, 9, <fc). Figure 12 illustrates how the physical obstacle 
has a new "shadow" or umbral region on y 2 for every 6 £ 0. A convenient 

method of dealing with this situation, which is complicated by the fact that the 
integration for E <x', y\ z') is actually taken on the surface y 2 and the shadow 
contours may be quite general in any practical case, is to introduce a reference 
plane (II K ). 1 

For a given observer position (°°, ( 9 , y ) , the area to be deleted from the integra- 
tion (0 2 ) is identified on y 2 by its boundary contour (C 2 ). These points can be 
projected into the reference plane (1^) to form 0 2 . Where a large number of 
obscurations are present this method is especially convenient, since all deletions 
can be accomplished by means of a Cartesian net on II K . The true area deleted 
on y 2 is related very simply to the projected area in II via 1 

AS = sec y* A A 

(increment of area on y 2 ) (increment of area on II K ) 


1 Reference 6, p. 208. 
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where the angle y* is obtained from the dot product between the unit normal (n 2 ) 
to the surface y 2 and the basis vector ft. 

cos y* - n 2 • k 


It is also possible to employ curve-fitting to the contour C 2 ' on II K in the deter- 
mination of AA. 

The preceding ideas were used in the present report to estimate the effects of 
main reflector seconciary obscuration by the hyperboloidal subreflector, the flat 
mounting plate, and the 14-inch diameter trusses. Primary obscuration was not 
included at this stage of the development. Multiple scattering was also ignored. 
Methods of estimating current distributions on support trusses, and computing 
the scattered fields of members whose dimensions in terms of a wavelengths are 
sometimes small, have been considered, but not implemented. Since the directive 
gain for the antenna was very high, and the beamwidth very small, all of the side- 
lobe structure of interest was observed for Q° < & < 0.50° at '15.3 GHz and 0° £ 

6 < 0.80° at 7.3 GHz. The obscurations or "shadows" due to subsystem and sup- 
port struts were therefore not updated with observer angle. 
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APPENDIX B 


A COMPARISON OF MEASURED DATA WITH ANALYTICAL RESULTS 


At the conclusion of the series of 15.3 GHz computations contained in the body of 
this report it was determined that only 7.3 GHz radiation-pattern data would be 
available for this analysis. Only the ''complete” configuration which included 
mounting plate , subsystem obscuration, and strut blockage effects was recom- 
puted at 7.3 GHz. The results were compared with measured data, 0 = 0° cut. 
obtained via DOD Satellite Number 16. Figure 13 shows the measured data, pro- 
vided by the U.S. Naval Research Laboratory. An undetermined amount of noise 
can be seen for signal levels approximately 20 db below the beam maximum. 

Since reduction of the simulation frequency from 15.3 GHz to 7.3 GHz eased 
the core and cpu requirements significantly, a valuable re-evaluation of the tech- 
niques used in the body of the report could be made. Two obvious areas of ex- 
ploration were considered. First, a study of the effects of interpolation of the 
backscattered fields from the subsystem composed of a hyperboloidal cap and a 
flat mounting plate was made. The theta ( 6 ) increment used in Figures 4 through 
7 was Ad =0.2 degrees implying approximately 40 data points in the domain of 
interest (100° £ 6 £ 180°). For a sampling criteria of LI 2 = 3.0, however, at a 
wavelength of 0.06431’ (15.3 GHz) the amount of data required along a 30-foot 
parabolic arc is approximately 300 points. Ordinarily this poses no special 
problem when the subsystem behaves as a virtual point source. Introduction of 
the flat annu lar mounting plate perturbed the spherical wavefront, see Figure 7 , 
therefore the required interpolation of field data is no longer empty. 

Second, a study of the assumption that the backscattered data from the hyper- 
boloidal cap and flat annular mounting plate could be computed. on the arc R = F, 
and then expanded as a spherical wave up to the paraboloidal surface, was under- 
taken. As mentioned previously, there is some justification for attempting to 
observe the subsystem fields on a circular arc. The spherical wavefront, or 
slight departures from it, can be displayed. A decision can then be made relative 
to the subsystem sampling (LIj), Furthermore, a parallel experimental effort 
would logically utilize a circular arc such as R = F to probe the subsystem fields. 
If ray optics is used to identify the virtual origin of the subsystem fields, see 
Figure 14, the conclusion is that the hyperboloidal cap has a unique phase center 
situated at the focal point F, but the flat annular mounting plate has a unique 
phase center a distance Z = FF* from the focal point. These phase centers 
have, in general, different amplitude weighting -factors and initial phases asso- 
ciated with them. 
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In view of the preceding, the expansion of the subsystem fields according to a 
spherical law, with phase center at focal-point F, is difficult to justify. Since 
one does not a priori know the location and weight of the various phase centers 
and expansion about multiple phase centers would be tedious, requiring eventual 
superposition at the paraboloid, it seemed expedient to compute the subsystem 
backscattered fields directly on the parabolic locus of the main reflector via 
diffraction theory. A parallel experimental effort should also proceed along the 
parabolic law which is to be illuminated, an undesirable complication which 
would require a separate measurement dependent on each paraboloidal F/D ratio 
to be employed. 

In order to identify the variations introduced by the alternative methods of com- 
puting the backscattered phase and amplitude of the fields from the compound 
subsystem, the 7.3 GHz computations were first done with interpolation on a 
spherical arc (R = F). The amount of data in the domain of interest (100° £ 8 < 
180°) was approximately 40 points. The amount of data required on the 30-foot 
parabolic arc for ^ = 0.13' (7.3 GHz) is approximately 100 with LI 2 set at 3.0 
for Kirchhoff integration. It can be seen that the interpolation is still significant 
but less critical than at 15.3 GHz. See Figure 15. Four vertical lines on Fig- 
ure 15 represent the geometrical bounds of the problem. The two outer bounds 
are the edge of the main reflector and the obscuration due to the subreflector, 
while the inner bounds show the extremes of the rays originating at the mounting 
plate under a ray-optics analysis. Figure 16 presents the far-field pattern for 
the complete Cassegrain system. This result exhibits a monotonic descent in the 
sidelobe structure which is not seen in the measured data. Due to the presence 
of noise in the measured pattern, the absolute signal level beyond the first side- 
lobe is open to question. 

At this point the 7.3 GHz calculations were repeated, but without interpolation of 
the baekscattered-field data. A spherical locus (R = F) was retained for these 
computations. Figure 17 shows the subsystem pattern by this approach. Each 
of 64 field values was obtained at a spherical angle ($) pre-computed by means 
of a separate driver program to correspond to the (a, £) net LI 2 = 3.0 criterion 
on the main reflector arc. Figure 17 should now be compared with Figure 15'. 

The domain of interest was contracted to 100° < 9 < 166° to conserve cpu time 
since the values outside of the extreme geometric bounds do not enter into the 
Kirchhoff integration. Significant differences appear between Figure 17 and 
Figure 15 due to avoidance of the interpolation process which was previously 
routinely carried out using E 0 and 4>z 0 analogue subsystem backscattered data 
from the GSFC program. 

The far-field pattern corresponding to Figure 17 is Figure 18, and should be 
compared with Figure 16. Evidently the avoidance of interpolation has improved 
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the quality of the result since the sidelobe envelope is no longer monotonically 
decreasing, but reproduces the prominent lobe whose peak lies in the vicinity of 
6 - 0.5 degrees in the measured pattern. It is unfortunate that the measured 
pattern of Figure 13 contained enough noise to mask nulls and raise sidelobe 
levels at the -25 db to -30 db level relative to the beam maximum. 

Next a parabolic locus, the arc of the main reflector, was used to compute 
the backscattered subsystem fields at 7.3 GHz. Interpolation was avoided as in 
the preceding case. Figure 19 shows the field amplitude arriving at the surface 
of the paraboloid prior to conversion to charge distribution and electric sheet 
current. The associated phase is not shown since phase varies rapidly on the 
parabolic arc and order is difficult to recognize with the inclusion of the effects 
of the flat annular plate. The corresponding far field pattern is Figure 20 and 
shows considerable departure from Figure 18 and Figure 16. Most interesting, 
perhaps, is the fact that Figure 20, the best estimate of the far-field pattern ob- 
tained with the diffraction program at 7.3 GHz, reproduces a first sidelobe peak 
at 6 - 0.24 degree and a second sidelobe peak at 6 ~ 0.49 degree as in the 
measured result, Figure 13. The level of the first sidelobe via the program is 
-20.5 db, and is in exact agreement with measured results. The level of the 
second sidelobe was computed as -27 db, and a meaningful quantitative com- 
parison at this level cannot be made since noise has spoiled the measured pat- 
tern beyond 0.30 degree as evidenced by obliteration of nulls beyond this value 
of theta. It is noted that the computed result, Figure 20, predicts substantial 
null-filling around d = 0.375’ degree due to the combined effects of the annular 
plate subsystem obscuration and strut blockage. 

A review of the three approaches used to obtain the 7.3 GHz diffraction pattern 
forces the conclusion that the parabolic locus approach without interpolation of 
the backscattered subsystem fields is the better approach. This is particularly 
true when there is a large departure from the virtual point-source that emanates 
an ideal set of spherical wavefronts. The favored approach is also the most de- 
manding in terms of cpu time and core requirements. A comparison of computed 
data with several USNRL patterns obtained by satellite at 7.3 GHz, although in 
the presence of noise interference, indicated that the spherical -wave expansion 
assumption and subsystem data interpolation were responsible for significant 
inaccuracy beyond the main beam and first sidelobe computation. The evidence 
appeared convincing enough to warrant a re-examination of the 15.3 GHz problem 
for the "complete" configuration. 

Since Figure 7 and Figure 11 in the body of this report present the backscattered 
subsystem pattern on a spherical locus using interpolation and the corresponding 
far-field pattern, respectively, the interpolation was avoided in a subsequent 
study. This produced Figure 21 for the subsystem and Figure 22 for the far-field 
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patterns. It is noted that although the hyperboioidal cap is confocal with the 
paraboloid and all rays incident upon it emerge phase-coherent in the paraboloid 
aperture-plane, theeffect of the flat annular plate, situated in the focal -plane,* 
is frequency dependent. The influence of the flat plate can once again be seen by 
noting the effect on phase and amplitude of the backscattered fields of Figure 21 
between the innermost geometric bounds of the flat plate. A comparison between 
Figure 7 and Figure 21, as well as between Figure 11 and Figure 22, shows the 
variations introduced by interpolation of the subsystem data. The principal 
features are very similar although details are seen to vary in the far-field 
patterns. 

Finally the parabolic locus was employed, without subsystem data interpolation . 
'Figure 23 and Figure 24 present the computed subsystem and far-field patterns, 
respectively. Satellite data were not available for comparison. This pattern, 
which represents the culmination of the entire effort, exhibits a first sidelobe 
level of -20 db at Q - 0.115 degree, obliteration or filling of the second null, and 
a second sidelobe level of -28 db at Q - 0.230 degree. The pattern at 15.3 GHz 
is characterized by featux'es very similar to those at 7.3 GHz, where the first 
and second sidelobe levels were -20 db and -27 db, and where the second null 
was also obliterated. Since measured data of high quality were not available at 
15.3 GHz it is necessary to appeal to the comparison between measured and 
computed results at 7.3 GHz and rely entirely on computed or simulated results 
at the higher frequency. 


*This would not be true if the plate were midway between the feed point F* and the focus F of 
the paraboloid as this would constitute a degenerate pair of coincident hyperboloids. 
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Figure 1. NRL Microwave Space Research Facility (Waldorf, Md.) 
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Figure 2. Far-field Radiation Pattern (f = 1.53 GHz) 
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Figure 4, Hyperboloid Radiation Pattern at R = F (JA 1 = 1,25) 
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Figure 6. Hyperboloid Radiation Pattern at R = F (LI t ■- 0.300) 
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Figure 9. Far-field Pattern Showing Hyperboloid Blockage Effect 
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Figure 10. Far-field Pattern Showing Hyperboloid-plate Diffraction 

and Blockage Effects 
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Figure 11. Far-field Pattern Showing Hyperboloid-plate Diffraction 
and Blockage Plus Truss-blockage Effects 
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Figure 13. 7.3 GHz Radiation Pattern (Measured) 
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Figure 14. Reflector-System Geometry 
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Figure 15. Hyperboloid-plate Radiation Pattern at R = F (LIj = 0.3) Using Interpolation 
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Figure 17. Hyperboloid-plate Radiation Pattern at R = F (LIj = 0.3) Without Interpolation 
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Figure 18. Far-field Pattern (f = 7.3 GHz) 
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Figure 21. Hyperboloid-plate Radiation Pattern at R = F (LI t = 0.30) Without Interpolation 
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Figure 22. Far-field Pattern (f = 15,3 GHz) 
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Figure 24. Far-field Pattern (f = 15.3 GHz) 
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